Fano resonance resulting from a tunable interaction between molecular vibrational 
modes and a double-continuum of a plasmonic metamolecule 
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Coupling between tuneable broadband modes of an array of plasmonic metamolecules and a 
vibrational mode of carbonyl bond of poly (methyl methacrylate) is shown experimentally to produce 
a Fano resonance, which can be tuned in situ by varying the polarization of incident light. The 
interaction between the plasmon modes and the molecular resonance is investigated using both 
rigorous electromagnetic calculations and a quantum mechanical model describing the quantum 
interference between a discrete state and two continua. The predictions of the quantum mechanical 
model are in good agreement with the experimental data and provide an intuitive interpretation, at 
the quantum level, of the plasmon-molecule coupling. 
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The spectral control of optical absorption in metallic 
nanostructures has generated rapidly growing interest, 
both as a means for studying light-matter interaction at 
the deep-subwavelength scale and to develop new func- 
tional nanomaterials. Plasmonic materials provide an 
ideal platform for achieving enhanced optical interaction 
at the nanoscale as their primary building blocks consist 
of highly resonant particles, which strongly confine and 
enhance the optical field [TH1] . These properties are com- 
monly employed to increase the optical coupling between 
plasmonic structures and an optically active medium, 
enabling many applications, including surface-enhanced 
Raman spectroscopy [5 , 6 , surface-enhanced infrared ab- 
sorption spectroscopy [7J [8], enhancement of the excita- 
tion rate and fluorescence of quantum dots (QDs) and 
molecules [9HTT]. chemical sensing [I2j [13], and biode- 
tection [14j [15] . One effective approach for tailoring the 
shape of plasmonic resonances is the Fano interference 
between a narrow (discrete) resonance and a broad (con- 
tinuous) one [T6HT9] . 

A defining feature of Fano resonances is their pro- 
nounced spectral asymmetry. This can be exploited in 
applications that require high sensitivity to changes in 
the environment. Fano resonances in plasmonic systems 
can also provide physical insights into light-matter inter- 
action at the nanoscale, as they are primarily the result 
of short-range, near-field interactions. Although Fano 
resonances in plasmonic systems have been chiefly inves- 
tigated in a classical context, they can also be observed in 
quantum systems consisting of QDs, atoms, or molecules 
coupled to metallic nanoparticles [14 , 20 - 22 . In this case 
the Fano resonance results from the interaction between 
a quantum discrete state and a classical continuum (lo- 
calized or extended) plasmon mode. In many applica- 
tions it is desirable to tune the strength of the inter- 
action between the discrete and continuous state, and 
consequently dynamically change the shape of Fano res- 



onances. However, since the materials and geometrical 
parameters of a plasmonic structure are difficult to con- 
trol at the nanoscale, a convenient way to achieve optical 
tunability has not yet been attained. 

In this Letter we demonstrate experimentally and the- 
oretically an efficient approach to achieve in situ tun- 
ing of the strength of the interaction between optically 
active vibrational modes of molecules - here the car- 
bonyl bond of poly(methyl methacrylate) (PMMA) - and 
the localized surface plasmon (LSP) modes of plasmonic 
metamolecules consisting of asymmetric cross apertures 
in a thin metallic film (see Fig. [I]). These nanostruc- 
tures have been studied extensively [23j [24] , one of their 
defining properties being that, analogous to anisotropic 
molecules, their polarizability strongly depends on the 
polarization of the incident wave. These plasmonic meta- 
molecules can therefore be tailored to different molecular 
resonances not only because the resonance frequencies of 
LSPs are strongly dependent upon their size and shape 
but, more importantly, because of their intrinsic optical 
anisotropy. 

Sample fabrication. — A CaF2 substrate of thickness 
1 mm was partially coated with a 5 nm chromium adhe- 
sion layer and a 35 nm gold film by thermal evaporation. 
The array shown in Figs. [TJa) and [ljb) was patterned 
using electron-beam lithography and argon-ion milling. 
The aperture geometry, and the orientation of the elec- 
tric field, E, of incident polarized light relative to the 
apertures is shown in Fig. [ljc). The apertures have 
arm lengths L x - 1900 nm and L y - 1340 nm and arm 
widths G x = 580 nm and G y = 360 nm. The surface was 
then spin-coated with a 85 nm PMMA film, filling the 
apertures. Figure [TJd) shows the final structure in cross 
section. 

Optical characterization. — The transmission and re- 
flection spectra of the array were measured using Fourier 
transform infrared spectroscopy (FTIR) microscopy. The 
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FIG. 1. (color online), (a) Scanning-electron micrograph 
showing part of the cruciform aperture array before PMMA 
coating, (b) magnified detail of four apertures, (c) aperture 
geometry, showing the in-plane electric field polarization an- 
gle, 0, (d) schematic cross section of the structure following 
PMMA coating. 



source is a broad-band mid-IR globar which is focused 
onto the sample. The light is collimated and passed 
through an aperture which selects radiation being trans- 
mitted or reflected from a circular sample area 79 /xm in 
diameter. The signal was detected using a HgCdTe detec- 
tor. The reflection spectra R and transmission spectra 
T of the samples were normalized to spectra obtained 
from uncoated areas of gold film and CaF 2 substrate, re- 
spectively. The incident radiation was linearly polarized 
using a wire grid polarizer with an extinction ratio of 
168 : 1 at 5 jam. 

Figures |5|a) and|2|c) show, respectively, the measured 
transmission and reflection spectra of the aperture ar- 
ray coated with PMMA. In the transmission spectra two 
global maxima are present; one, at 4.2 jam occurs when 
the polarization angle 6 = and another, at 5.5 /im, when 
= 90°. These global maxima are due to the excitation 
of two LSP modes, which are supported by the aper- 
tures [23j |24] . In the reflection spectra the excitation of 
the plasmon modes produces global minima. The spec- 
tral signature of the PMMA carbonyl bond is located 
at 5.79 /im. The shape of the carbonyl resonance peak 
changes with showing a more pronounced, asymmetric 
Fano-like line shape as approaches zero. The interac- 
tion of the two broad plasmon resonances with the nar- 
row molecular resonance forms the Fano resonance. The 
inset to Fig. [2|a) also shows that, for small values of 0, 
at wavelengths somewhat longer than the carbonyl bond 
resonant wavelength, the transmission can be larger than 
the baseline value in spite of the fact that PMMA is a 
lossy dielectric. This is an example of absorption-induced 



FIG. 2. (color online), (a) Transmission and (c) reflection of 
the structure measured using FTIR; (b) transmission and (d) 
reflection simulated using the RCWA method. The incident 
polarization angle varies from zero (brown) to = 90° (pur- 
ple). The resonance of the PMMA carbonyl bond is visible 
at A = 5.79 /xm. The inset to (a) shows the measured spec- 
trum (solid line) and the fitted base-line (broken line) at the 
carbonyl bond resonance for = 0°. 



transparency [25]. Furthermore, as shown in Fig. [2|c) 
the reflectivity of the aperture array decreases with re- 
spect to the baseline for small whereas it increases for 
large (9, by as much as -100% for = 90°. This phe- 
nomenon we refer to as absorption-induced reflectivity. 

Numerical simulations. — To model the optical re- 
sponse of our PMMA-coated device we used commer- 
cially available software [26] based on the rigorous cou- 
pled wave analysis (RCWA). We used a constant refrac- 
tive index for the CaF2 substrate, a Lorentz-Drude model 
for the Au and Cr layers [27 , and a complex refractive 
index for the PMMA layer (found by fitting the numer- 
ically determined transmission to an experimental mea- 
surement of the transmission of a 85 nm thick PMMA 
film). Figures |2|b) and|2|d) show the transmission and 
reflection spectra, respectively, found by numerical sim- 
ulations; there is good agreement with the experimental 
data. The simulations also allow spatial field profiles to 
be calculated; Fig. [3] shows the spatial profile of the field 
enhancement corresponding to the two LSPs. 

Theoretical model for molecule-plasmon metamolecule 
interaction. — In what follows we introduce a simple quan- 
tum mechanical model describing the interaction between 
the molecular vibrational mode and the LSP resonances. 
Figure [3] shows the interacting quantum mechanical sys- 
tem schematically. 

We assume that the interacting system is described by 
a Hamiltonian, H = Ho + V a + Vb, where Ho = H m +H a +Hi) 
is the Hamiltonian of the noninteracting molecule (ff m ) 
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FIG. 3. (color online). System schematic showing the ground 
state \g), the two LSP continua, \E;a) and \E;b), and the 
discrete molecular resonance |e). Simulated spatial profiles 
(taken in the x - y plane at the center of the aperture) show 
large field enhancement for # = 0atA = 4.2 /xm and = 90° 
at A = 5.5 /im. Aperture boundaries are indicated by dashed 
lines. The atomic structure of the PMMA monomer is shown 
with the carbonyl bond in pink; the resonance of this bond 
corresponds to the excited state |e). 



and the two LSPs (H a and Hb) and V a and V b describe 
the interaction between the molecule and the two LSPs. 
The molecule is modeled as a two- level system with \g) 
and |e) the ground and excited states, respectively. We 
assume that the eigenstates corresponding to the energy 
E, \E\a) (\E\b)), of H a (Hb) form two continua. Note 
that the two plasmon modes are orthogonal and thus the 
continua are orthogonal, too: {E; a\E f ; b) = S a bS(E - E r ). 

Under these conditions we can write Hq as 



H = E g \g)(g\+E e \e}(e\+ £ f E\E;i)(E;i\dE, 

i=a,b 



(i) 



where E g and E e are the energies of the ground and ex- 
cited state of the molecule, respectively. The interac- 
tion between the molecule and the two LSPs (continua) 
can be expressed as Vi = f Vi^\E] i)(e\dE + h. c, i = a, 6, 
where V^e = (E]i\Vi\e),i ~ a > ^ are the matrix elements 
of Vi and are proportional to the local electric field of the 
LSP and h. c. means Hermitian conjugate. Introducing 
the density of states of the two continua, p aj b(E), these 
matrix elements are written as Vi,e - viyj pi(E),i = a, b. 
We assume that the spectra of the plasmon A (B) is de- 
scribed by a Lorentzian centered at the plasmon energy 
E a (E b ) with width T a (T 6 ), 



Pi(E) 



1 



: a, b. 



(2) 



Following the approach of Fano [16 , the total Hamil- 
tonian H can be orthogonalized in a new basis of 
eigenstates which form two orthogonal continua, i.e., 
(E;l\H\E';l) = (E;2\H\E f ;2) = E5(E - E'), where 
(E; l\E f ; 2) = 0. These eigenstates can be written as [16] 



\E;l) 



sin A 
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^)-cosA|^); \E;2) = \^ E ), (3) 



where | $>e) is the "dressed" excited state |e) modified by 
the interaction with the two continua, 



;> = |e>+ £ f 



dE' 



yfin(&j(E - E>) 



\E';i), (4) 



(here f means the Cauchy principal value of the integral) 
and \i(j e ) and \i(je) are two continuum states defined as: 



= y/*h{E) (V a , E \E; a) + V b , E \E; b)) , (5a) 
$ E ) = ^H(E) (V b * E \E; a) + V a *jE; b)) . (5b) 

In Eq. ([3| A determines the phase shift induced by quan- 
tum interference between the discrete state |e) and the 
two continua, and 7 describes the width of the modified 
excited state. They are given by: 



7(£) = tt(K,£| 2 + |H,£| 2 ) 



A(E) 



- arctan - 



1{E) 



E-E e -X(E)' 
where the energy shift of the excited state is 

l(E') 



E(E) = f - 

J 7T 



(E-E') 



dE'. 



(6a) 
(6b) 

(7) 



The interaction between the coupled molecule- 
metamolecule system and a plane wave normally incident 
onto the array can be described by a Hamiltonian T = 
f m + T a + f 6 , with f = teo\e)(0\ + £*=a,6 / T itE \E] i){0\dE + 
h. a. Here, |0) is the ground state of the interacting sys- 
tem, in which the molecule is in the ground state and 
no plasmons are excited. The matrix elements of T are: 
t e o = <e|T m [0), T a , E = ( E-a\f a \0) = p a ^p a (E)smO, and 
Tb,E = (E;b\f b \0) = pb\/ 'pb(E)cos6, where p ajb quantify 
the coupling strength between the incoming plane wave 
and the two LSPs and we have used the fact that the 
plasmon A (B) is excited by a ^/(x)-polarized wave. 

Using Fermi's golden rule, the rate at which the 
system absorbs photons of frequency uj is W(oj) = 
^Y.i=i,2l dE\(E]i\f\())\ 2 5(huj - E). Inserting Eqs. 
in this equation one obtains the absorption of the system, 

A(w) = Afco [\& h jT\0)\ 2 F(q, e) + |<VUT |0)| 2 ] , (8) 

where J\f = 27rn s /io is a constant, with n s and Iq being 
the number of molecules per area and the intensity of the 
incoming beam, respectively, and F(q, e) = (<2+e) 2 /(e 2 +l) 
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is the Fano function. The asymmetry parameter, q, and 
the reduced energy, e, are given by 



0.5 



<$UT|0) 



V7T7(M<^Jr|o)' 



huo - E e - E(ftcj) 



(9) 



The absorption spectrum in Eq. ([8| shows that the 
coupling of the molecular vibrational mode with the two 
continua of the plasmonic metamolecule leads to a res- 
onant response of the optical system around the energy 
of the vibrational mode, huj - E e , which is superimposed 
on a broad background. Hence, unlike the case of only 
one continuum, the absorption does not cancel at the 
frequency for which F(q,e) = 0, i.e., when q = -e. 

The interaction between the incoming light and the 
plasmons can be tuned by varying the polarization an- 
gle; this allows one to modify the mixture of the two LSPs 
and, consequently, one can easily tune the strength of the 
interaction between the vibrational mode and the plas- 
monic metamolecule. In order to provide a more quanti- 
tative description we assume that the normalized densi- 
ties of state of the two LSPs are given by Eq. Then, 
the matrix elements in Eq. ^ and the parameters of the 
Fano function can be calculated explicitly. The absorp- 
tion can be cast as 



A(oj) 



AT* 



[ifJLaVaPa sin # + ^bPb cos#| 2 F(<?, e)n 



PaPbl^aVb sin + p h v a cos 0\ 2 ] , (10) 
where the Fano parameters are 
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-^{pa\v a \ 



+ Pb\Vb 



(11a) 
2 huj - E b 
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Comparison between theory and experiment. — To fit 
the model to the experimental system the absorption A of 
the structure was extracted from the measured transmis- 
sion and reflection data using the relation A = 1 - T - R, 
as shown in Fig. Qa). From this the parameters E a ^ and 
r a> 5 can be directly extracted; E a and T a from the 6 = 
data and E b and from the = 90° data. E e and T m 
(the carbonyl bond parameters) are taken from trans- 
mission measurements of a 100 nm PMMA film. The 
parameters v a ^ : /i a? &, and t e o are then extracted for each 
value of 6 using an unconstrained least squares method. 

Figure [4jb) shows the predictions of the model follow- 
ing fitting. There is good agreement between the mea- 
sured spectra and those predicted by the model at all val- 
ues of 0. In particular, all the main features of the Fano 
resonance are reproduced by our theoretical model, e.g. 
the increase of the optical absorption as is varied from 
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FIG. 4. (color online), (a) Absorption of the structure, from 
FTIR measurements of transmissio n a nd reflection, (b) Best 
fit of experimental data using Eq. (10). 



zero to 90°. Our model also predicts that t e o « \v a ^\ 2 for 
all angles #, which means that, as expected, the molecules 
interact primarily with the optical near-field of the SPPs. 
This proves that our model can be used to gain insights 
into the interaction between quantum and plasmonic sys- 
tems as it allows one to retrieve the values of basic quan- 
tities, e.g. v a ^, /i tt5 6, and t e o, which characterize these 
interactions at the nano-scale. 

In conclusion, we have demonstrated that plasmonic 
structures allow in situ tuning of the optical coupling 
between a quantum and a classical system simply by se- 
lecting the polarization of incident light. This tuning 
produces a dramatic change in the spectral line shape 
about the optical resonance of the quantum system, a 
consequence of the Fano resonance produced by the cou- 
pling of broad and narrow resonances. The system has 
been described using an intuitive quantum mechanical 
model, which allows a deeper understanding of the un- 
derlying physics. Our study has the potential to foster 
exciting new developments in sub-wavelength optics and 
nanophotonics and in the emerging field of quantum plas- 
monics. Equally important, the system investigated here 
has the potential to improve the range and versatility of 
a series of nanodevices, including chemical and biological 
plasmonic sensors. 
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